Determination of the exchange energies in Li 2 VOSi04 from a high-temperature series 
analysis of the square lattice J\-Ji Heisenberg model. 
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We present a high-temperature expansion (HTE) of the magnetic susceptibility and specific heat 
data of Melzi et al. on Li 2 VOSi0 4 [Phys. Rev. B 64, 024409 (2001)]. The data are very well 
reproduced by the J1-J2 Heisenberg model on the square lattice with exchange energies Ji = 1.25 ± 
0.5 K and J2 = 5.95 ±0.2 K. The maximum of the specific heat C^^Tmax) is obtained as a function 
J2/J1 from an improved method based on HTE. 



PACS numbers: 75.10.Jm 75.50.Eo 

I. INTRODUCTION 

The vanadium oxide Li2VOSi04 is a quasi two- 
dimensional magnet which is well described by the spin- 
i Heisenberg model on the square lattice with first (Ji) 
and second (J2) neighbors interactions : 

h = ^ s l -s J + j 2 3 ■ ^ w 

<i,j> <<i.j>> 

In their NMR experiments Melzi et alS^ found evidence 
for a Q — (it, 0) or (0, it) ordering in this system, as 
expected in the classical spin model for J 2 > 0.5 Ji (order 
from disorder) and also as expected in the spin-i case for 
</2^0.6Ji£ From the high temperature behavior of x(T) 
the Curie- Weiss temperature 6 = J\ + J 2 was estimated^ 
to be ~ 8.2 ± 1 K. Combining this information with 
the position of the maximum of the specific heat of the 
system (and comparing it to exact diagonalization data 
for 16 spins), Melzi et alA estimated the ratio of exchange 
energies to be J2/J1 — 1.1. 

In two recent papers, Rosncr et al£& considered 
the determination of J± and J 2 in Li2VOSi04. They 
performed band-structure calculations (local-density ap- 
proximation) and found J\ + J2 — 9.5 ± 1.5 K and 
J2/J1 — 12. They also computed the high-temperature 
expansion (HTE) of the magnetic susceptibility and of 
the specific heat for the Hamiltonian of Eq. ^ The x(T) 
and C V (T) obtained from these series at J2/J1 — 12 are 
in reasonable agreement with the experimental data for 
\(T) and C V (T). Rosner et al£ did not use the HTE to fit 
the data in the high temperature region (T > 20 K) be- 
cause they could not find a Curie law behavior. Melzi et 
al. indeed had to subtract a temperature-independent 
constant to recover a 1/T behavior at high temperature. 



This term is due to the Van-Vleck paramagnetism (spin- 
orbit coupling effect) and its magnitude is typical for a 
V A+ ion in a pyramidal environment^ It is of course cru- 
cial to exploit the high temperature region of the ex- 
perimental data in order to take full advantage of the 
HTE. In this work we show that the HTE for the sus- 
ceptibility allows us to determine the exchange energies 
in Li2VOSi04. We show that this method is an un- 
biased way of determining the microscopic parameters 
of this model. As a result we find a ratio J2/J1 — 5 
which is significantly different from the value 12 predicted 
previously^ We checked this result by a calculation of 
the specific heat. For this purpose we employed a recent 
technique^ which extends the convergence of the series 
to low temperatures by using two constraints related to 
the total entropy of the system and to its ground-state 
energy. 



II. HIGH TEMPERATURE SERIES AND 
EXTRAPOLATION 

The series for x(T) and C V (T) are computed in the 
standard way by a cluster expansion method to order 
1/T . Each coefficient is a rational number and is com- 
puted exactly (series available upon request). These 
series were computed independently in Ref. |6j to order 
1/T 10 and our results agree with their. 



III. FIT OF THE SUSCEPTIBILITY 

In the first step we fit the susceptibility data of Melzi 
et alrk by the following procedure: 
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• Temperature range of the fit. The experimen- 
tal data are fitted over a fixed temperature range 
[T min ,T max ] where T max = 300 K is the highest 
available temperature available in Ref. |2j and T m ; n 
is determined self-consistently so that the suscepti- 
bilities calculated from the HTE are well converged 
down to T m ; n when J% and J 2 are close to the op- 
timal parameters. A rather high T m ; n = 8 K is 
chosen for the first global scan of Fig. ^ this in- 
sures that the calculated x(T) is well converged in 
the whole interval whatever (J l7 J 2 ). Once the rel- 
evant region in J1-J2 space is identified, T min is 
lowered down to T mnl = 5.5 K, a temperature still 
above that where the various approximants approx- 
imants (Pade) start to differ from each other. 



• Scan over Ji and J 2 . For each set of couplings, 
the HTE for *(T) (including 1/T 11 order) is ex- 
trapolated by all possible Pade approximants. The 
rational fractions which have (spurious) poles in the 
interval [T m i„,T max ] are discarded. 



• Measure of the error. The experimental data 
X exp are compared to xw + Co X th (T) where 
is Xw is a temperature-independent Van-Vleck 
term and Co the Curie constant. For each 
value of (Ji,Ja) and for each Pade the pa- 
rameters xw an d Co are determined to insure 
the best possible fit. The quality of the fit 
is measured in a standard way through S 2 = 



1 6 ! 
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FIG. 1: Quality of the susceptibility fit as a function of Ji 
and J2. The biggest circles correspond best fits. Two regions 
around J2 ~ 6 K and Ji ~ 1 K (resp. J2 ~ 1 K and Ji ~ 6 K) 
appear to be compatible with the data on Li2VOSi04. The 
radius of the circles are proportional tor = 8E^ in — E 2 where 
E 2 is the mean square difference between the theoretical curve 
and experimental data and E^in corresponds to the optimal 
(smallest value of E 2 ). Parameters which give r < are not 
represented. Dashed circles : J Cv — 14.5±1.5 K 2 . The dashed 
line corresponds to a Curie- Weiss temperature Ji + J 2 = 7 K. 



IV. SPECIFIC HEAT 

At high temperature, the specific heat of the J1-J2 
model behaves as: 



The figure ^ shows that two regions of the parameter 
space are compatible with the susceptibility of the J1-J2 
model. The points where E 2 is larger than eight times 
its smallest value (fit of poor quality) are not displayed. 
It turns out that the Pade approximants to the suscepti- 
bility converge down to T ~ 5.5 K in the two "patches" 
of Fig. ^ For this reason T m ; n can be lowered down to 
5.5 K to refine the scans. Although both regions give fits 
of comparable quality, the existence of (jr, 0) magnetic or- 
dering in Li 2 VOSi04 as well as the results of Ref. make 
the J\ > J 2 scenario extremely unlikely. For this reason 
we will only focus on the region Ji~1.25 K and J2—6 K 
in the following. The refined scan is shown in Fig. [5] The 
susceptibility obtained from the HTE is compared with 
the experimental results in Fig. |3| The Van-Vleck sus- 
ceptibility and Curie constant are parameters of the fit 
and the optimal values are in complete agreement with 
Ref. 0. In order to estimate the quality of the fit in the 
high-temperature region where the HTE is almost exact, 
we subtracted the 1/T Curie term from the theoretical 
curve as well as from the experimental data. Even when 
the leading Curie term is subtracted (bottom of Fig.[3J), 
the deviation of the Pade approximants from the exper- 
imental results remains within experimental error bars. 
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Although the specific heat data are not very accurate at 
high temperatures due to the subtraction of the phonon 
contribution, a plot of C V T 2 shows (see inset of Fig. 
that Jq ~14.5 ± 1.5 K 2 . This constrains the values of 
Ji and J 2 to lie between the two large circles in Fig. ^ 
which is compatible with the independent fit performed 
on the susceptibility. 

The specific heat is obtained from the HTE as de- 
scribed in Ref. and will not be reviewed here. In addi- 
tion to the series itself, this method requires the knowl- 
edge of three quantities: 

• The total entropy, that is ln(2) per site. This means 
f™C v {T)/TdT = Nk B ln(2). 

• Ground-state energy per site eo of the Hamiltonian. 
Since the energy of Eq^is zero at infinite temper- 
ature we have C v (T)dT = -Ne . 

• The low-temperature behavior of C V {T). It behaves 
as ~ T 2 (resp. ~ cxp(— A/T)) when the ground- 
state is Neel long-ranged ordered (resp. gapped). 
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FIG. 2: Quality of the susceptibility fit as a function of Ji and 
Ji. Same as Fig. Q When several concentric circles are visible, 
some Pade approximants are significantly different from the 
others at low temperature. The points fall in the rectangle 
defined by Ji + 3.25J 2 = 20.59 ±0.11 K and J2 — 3.25Ji = 
1.9 ± 1.8 K. 
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FIG. 3: Top: Experimental susceptibility (crosses) compared 
to the HTE calculation (full lines). The different Pade ap- 
proximants (degrees [u, d] indicated) are not distinguishable 
at the scale of the figure for T > 7 K. The best fit is obtained 
for: Curie constant Co = 0.333 emu/mol, Van-Vleck suscepti- 
bility Xvv = 4.014 x 10~ 4 emu/mol. Bottom: Deviation from 
the Curie law (same data as top panel). 



The ground-state energy of the first-neighbor Hcisen- 
berg antiferromagnet on the square lattice is known 
very accurately from quantum Monte Carlo simulations 8 : 
eo = — 0.6694Ji. This result also applies to the infinite- 
J2 limit of the J1-J2 model. In the frustrated case 
the ground-state energy has been computed by differ- 
ent numerical techniques, such as zero-temperature se- 
ries expansion^ From these result we constructed a sim- 



FIG. 4: Specific heat. Full lines: theoretical prediction for 
Ji = 1.25 K and J2 = 5.95 K. The different lines correspond 
to different Pades approximants. Their dispersion is a mea- 
sure of the error bars. Dotted line: experimental data. Inset: 
experimental specific heat multiplied by T 2 . We estimate that 
C v T 2 /Nk B -> 14.5 ± 1.5 K 2 (horizontal lines). The agree- 
ment between the calculated C v and the experimental data 
can easily be improved (data not shown) at high temperatures 
by changing slightly the T 3 contribution (phonons) that was 
subtracted from the raw data^ 

pie ansatz which interpolates between the pure-Ji and 
pure-J2 models : 



e = Q'i J\ + a%J% - J (an J\ + a^Jf + a \2J\-h) (3) 

with a, x = -0.3135, a 2 = -0.1207, a u = 0.1267, 
«22 = 0.3011 and ai2 = —0.3722. The comparison be- 
tween the different approximants for the specific heat 
(with C v ~ T 2 ) at Ji = 1.25 K and J 2 = 5.95 K and the 
experimental data is shown Fig.0J Unlike usual Pade ap- 
proximants directly constructed on the HTE of the spe- 
cific heat, our procedure^ leads to accurate results (with 
a relative error smaller than a few percents) down to zero 
temperature. The agreement with the experimental data 
is very good for T > 4.5 K (no adjustable parameter). 
On the other hand we checked that the flat maximum 
observed in the experiment between 3 K and 4 K can 
hardly be accounted for by the J\-Ji Heisenberg model, 
whatever the couplings. This is a serious indication of a 
structural distortion (above T c ~ 2.8 K corresponding to 
the (3D) magnetic ordering of the system), as suggested 
by Melzi et a& from the analysis of the NMR spectra and 
by Becca and MilaiS from the theoretical point of view. 

When one of the J's is much larger than the other, the 
specific heat and the magnetic susceptibility are found 
to be nearly symmetric under the exchange J\ «-> J2 (at 
least at no too low temperatures). Thus this analysis 
cannot strictly discriminate between the two patches of 
Fig. ^ Nevertheless J2 > J\ gives slightly better fits. 

C V (T) for arbitrary Jij J\- — The specific heat is com- 
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FIG. 5: Maximum of the specific heat and temperature of 
that maximum as functions of J2/J1. The different curves 
correspond to the different Pade approximants at order 1/T . 
The regions where the model is assumed to have a ~ T 2 (resp. 
~ e _A//T ) specific heat are indicated. 

puted by the method above in a large range of coupling 
J2/J1. The behavior of C v at low temperatures is as- 
sumed to be: i) C V ^T 2 for J2/</i^0.38 (corresponding 
to an ordered antiferromagnet) ii) C v ~ exp(— A/T) for 
0.38< J2I Ji < 0.6 (gapped region^) and iii) C v ~ T 2 for 
0.6<J2/Ji (ordered antiferromagnet). The ground-state 
energy is approximated by Eq.|21 The results are summa- 
rized in Fig. [3] where the the value G™ lax of the maximum 
of the specific heat and the temperature T max at which 
C v reaches its maximum are displayed. The dispersion of 
the different approximants is again an estimate of the er- 
ror. The discontinuities in the curves are due to the fact 



that the Pade approximants which develop zeros or poles 
in the physical energy interval cannot be considered^ 
The results for C™ ax and T max represent a significant im- 
provement over the estimate made in Ref.0. The shape of 
the specific heat appears to be relatively independent of 
J2/J1 when J2 ^ 2Ji but some interesting structure ap- 
pear in the strongly frustrated region around J 2 — 0.5 Ji. 
The existence of a low C™ ax is indeed compensated by a 
small T max to conserve the total entropy. This shift of 
the entropy to lower temperatures and lowed energies is 
indeed expected close to quantum phase transitions and 
in frustrated systems in general. 



V. CONCLUSIONS 

We have computed the HTE for the uniform suscepti- 
bility of the J1-J2 Heisenberg model on the square lattice 
up to order 1/T 11 and obtained J\ = 1.25 ± 0.5 K and 
J2 = 5.95 ± 0.2 K by fitting the experimental data of 
L^VOSiO,}. The HTE for the specific heat to the same 
order and an improved method allowed us to extrapo- 
late C V (T) down to T = 0. A good agreement is found 
with the experimental data of L^VOSiC^for T > 4.5 K. 
Below that temperature the specific heat slightly differs 
from our results on J1-J2 Heisenberg model, which is 
consistent with the lattice distortion transition observed 
by Melzi et alA Chandra, Coleman and Larkinii pre- 
dicted an Ising-like finite temperature phase transition 
in the classical and quantum J1-J2 Heisenberg models. 
We could not find such a transition in the present cal- 
culations on the quantum model but some work is in 
progress to investigate thoroughly this question^ 
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